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^— ^ Abstract. We study the electromagnetic wave equation and the perturbed 

^— i massless Dirac equation on R< X R 3 : 

CN| ' utt - (V + iA(x)) 2 u + B(x)u = 0, iu t - Dm + V(aj)ii = 

5-H ■ where the potentials A(x), B(x), V (x) are assumed to be small but may be 

O-H rough. For both equations, we prove the expected time decay rate of the 

solution 

< ~II/IU 

7— H ■ where the norm can t>e expressed as the weighted L 2 norm of a few 

derivatives of the data /. 

c-j . 1- Introduction 

Dispersive properties of evolution equations play a crucial role in the study of 
nonlinear problems, and for this reason they have attracted a great deal of attention 
!— — in recent years. In particular, for the Schrodinger and the wave equation a well 
established theory exists, see [14] and [22]. On the other hand, in the variable 
coefficient case the theory is very far from complete. The simplest situtation is a 
| perturbation with a term of order zero; this is already very interesting from the 

■ physical point of view (electrostatic potential) . Several results are available for the 

equations 

id t u ~Au + V(x)u = 0, Uu + V{x)u = 0. 

I/-) , We cite among the others [8], [15], [16], [19], [32] and the recent survey [33] for Schro- 

dinger; and [5], [6], [10], [12], [13] for the wave equation. We must also mention the 
wave operator approach of Yajima (see [2], [39], [40], [41]) which permits to deal 
with the above equations in a unified way, although under nonoptimal assumptions 

Son the potential in dimensions 1 and 3. 
The next step in generality is a first order perturbation; from the physical point 
J> , of view this corresponds to a magnetic potential. In this case only a handful of 

results are available: Strichartz estimates for the 3D wave equation [11], provided 
the coefficients are small and in the Schwartz class; and smoothing estimates for 
the 3D Schrodinger and wave operators [37]. The most general case of variable 
coefficients has been studied in [17], [31] and [35], where local Strichartz estimates 
have been proved, in various degrees of complexity; see also [7]. 

In the present paper, our main focus will be on the three dimensional wave 
equation with an electromagnetic potential 

(1.1) u tt - (V + iA(x)) 2 u + B(x)u = 0, u:RxM 3 ^C, 

and the closely related massless Dirac system with a potential: 



(1.2) iu t -Vu + V(x)u = 0, 



n3 , ,p4 
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Here A : M 3 -> R 3 , B : R 3 -> R, V(a;) = V*(ar) is a 4x4 complex matrix on R 3 , 
and the symbol I? denotes the constant coefficient, elliptic, L 2 selfadjoint operator 

1 3 

3 = 1 

where the Dirac matrices ct\, a 2 , «3 have the following structure: 
(1.3) 
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We neglect the physical constants (i.e., we set c = H = 1), and we consider the zero 
mass case exclusively; the case of a positive mass, whose second order counterpart 
is the Klein-Gordon equation, has an additional term a±u with 



(1.4) 
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The relation between massless Dirac and wave equation is readily explained: indeed, 
the Dirac matrices satisfy the commutation rules 

aeatk + auai = 2<5 fc; / 4 

which imply immediately 

V 2 = -AI 4 , 

where I4 is the 4x4 identity matrix. Thus we have the fundamental relation 

(id t -V)(id t + V) = (A-d 2 t )h, 

which can be intepreted as follows: squaring the Dirac system produces a diagonal 
system of wave equations (or, conversely: taking the square root of a wave equation 
produces a Dirac system. According to the folklore, this was the route that lead 
Dirac to his equation). When a potential is present in the Dirac system, the above 
reduction produces an electromagnetic wave equation in a natural way. A discussion 
of this can be found e.g. in [23] (Volume 4, Chapter 4); see also section 6 below. 

Our goal here is to establish the decay rate of the spatial L°° norm of the solution, 
with minimal assumptions on the potentials. The expected decay rate is i -1 , both 
for the wave equation and the Dirac system. Indeed, known results for hyperbolic 
systems (for constant coefficients see e.g. [24], [25], and for C§° perturbations 
thereof see [20]) suggest a i^^ - decay rate in n space dimensions. 

Before stating our first result we introduce some basic notations. Under the 
assumptions of Theorem 1.1 below, the perturbed laplacian 

(1.5) H :=-{V + iA{x) f +B{x), 

where A(x) = (A^x), A 2 (x), A 3 (x)) : K 3 -> K 3 and B{x) : R 3 -> R, is a selfadjoint 
unbounded operator on R 3 ; the explicit standard construction is recalled in Section 
2. Spectral calculus allows us to define the operators tp(H) for any well behaved 
function ip( s )- 

In particular, consider a (non-homogeneous) Paley-Littlewood partition of unity 
on R 3 , defined as follows: fix a radial nonnegative function tp(r) € with i[)(r) — 
1 for r < 1, V(r) = for r > 2, define ^ (r) = V(2^ +2 r) - V>(2^ +1 r) for all j > 1, 
and tfio — ip. Then 1 = X)j>o *A? ^ s ^ ne required partition of unity on R 3 . The 
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operators (j>j(y/H) will be used in the following to define suitable norms associated 
to the operator H. We shall also use the notations 

(x) = (1 + M 2 ) 1/2 , (D) s .f = (1 - A) S/2 J = ^-\(0 s f) 
Our first result concerns the Cauchy problem for the wave equation perturbed 
with a small rough electromagnetic potential 

(1.6) u tt (t,x) - (V +iA(x)fu + B(x)u = 0, (i,i)£RxI 3 , 

(1.7) u(0,x) = 0, u t (0, x) = g(x). 
We can prove: 

Theorem 1.1. Assume the potentials A(x) £ K 3 , B(x) £ K satisfy 

3 



xl^Xllogl^H + 1)/=" Jl 1 laffl log + 1)0 ' 

/or some constant C > sufficiently small and some [3 > 1. TTiera any solution of 
the Cauchy problem (1.6), (1.7) satisfies the decay estimate 

(1.9) K^)l < y^2 2 -'||( a; )^ /2 ^(VH).g|| L2 , 

w/iere wp{x) := \x\(\ log |x|| + l)^ 3 . // m addition we assume that, for some e > 0, 

(1.10) (£>) 1+ %- e i 00 , (£>) e s e i°° 

i/ien u satisfies for any S > f/ie estimate 

(1.11) Ki,x)|<^||(.T) 3 / 2 + 5 .g|| H2+e . 

Remark 1.1. The norm appearing in (1.9) can be regarded as a distorted analogue 
of a standard Besov norm, generated by the operator Similar norms already 
appeared in [11] for magnetic potentials with coefficients in the Schwartz class; in 
that case, however, it was possible to prove the equivalence with standard Besov 
norms (see also [12], [13] for the analogous norms generated by — A+V(x), which are 
also equivalent to the nondistorted norms). Under the slightly stronger assumptions 
(1.10) on the coefficients, it is possible to prove an estimate like (1.11) expressed in 
terms of standard weighted Sobolev norms. 

Moreover, we remark that in our estimates we lose 2 derivatives; it is natural 
to conjecture that this is not optimal, and it should be possible to lose only one 
derivative as in the case of the free wave equation. 

Remark 1.2. As an essential step in the proof of Theorem 1.1, we need to establish 
the limiting absorption principle (LAP) for the operator H. This is obtained in 
Section 3 through several steps: starting from the "weak" LAP of [4] for the free 
resolvent, we first prove a strong version of the LAP for the free operator in the 
weighted spaces 

L 2 {w (j {x)dx), w fj (x) := |x|(|log|x|| + l)' 3 

and then we get the LAP for the perturbed operator. For the precise statements 
see Proposition 3.4. See also [37] for related results. 

Remark 1.3. When the initial data are of the form 

u(0,x)=f, u t (0,x)=0, 
Theorem 1.1 implies, by standard arguments, the estimate 

(1.12) \u(t,x)\ < j^&Mx^fniT/Hmv 

j>0 
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with an additional loss of one derivatives as expected. If in addition we assume 
that for some e > 

(1.13) (D) 2+c Aj g L°°, (D) 1+e B e L°° 
then also the simpler estimate 

(1.14) \u(t,x)\<j\\(x)^ 2+s f\\ H3+ ,. 
holds for all 8 > 0. 

Our second result concerns the perturbed Dirac system 

(1.15) iu t - Vu + V(x)u = 0, (t,i)elxl 3 , 

(1.16) u(0, a;) = /(*)• 

By explointing the above mentioned relation between the magnetic wave equation 
and the Dirac system, we can prove the following Theorem as a direct consequence 
of Theorem 1.1: 

Theorem 1.2. Assume the 4*4 complex valued matrix V(x) = V*(x) satisfies 

(!- 17 ) \V(*)\< , n wn f°i n , rwf . 1^)1 < 



N(x)(|iogN| + i)/3' 1 Wl - 1^(1^1x11 + iy> 

for some Co > small enough and some (3 > 1 . TTien the solution of the Cauchy 
problem (1.15), (1.16) satisfies the decay estimate 

(1.18) Mt,*)! < jJ2^\\(^T^^ + v)fh^ 

j>0 

where wp{x) = \x\(\ log |x|| + l) 13 . If in addition we assume that, for some e > 0, 

(1.19) (D) 2+e V e L°° , 
then u satisfies for any 8 > the estimate 

(1.20) \u(t,x)\<j\\(x) 3 / 2+5 f\\ H3+e . 

Since Theorem 1.2 is proved essentially by "squaring" the perturbed Dirac op- 
erator, a condition on the derivative DV is essential in order to apply Theorem 
1.1 to the resulting wave equation. On the other hand, we can study the Cauchy 
problem (1.15), (1.16) by a direct application of the spectral calculus for the self- 
adjoint operator T> + V(x); this alternative approach allows us to consider much 
rougher potentials V(x) (see (1.21)). The price to pay is an additional loss of one 
derivative, so that the total loss is 4 derivatives in our last result: 

Theorem 1.3. Assume the 4*4 complex valued matrix V(x) = V*(x) satisfies 

Co 

(L21) 1^)1 ^ | a; |l/2( a; )3/2(| log | a .|| + 1 )/3/2' 

for some Co > small enough and some f3 > 1. Then the solution of the Cauchy 
problem (1.15), (1.16) satisfies for any e > the decay estimate 

(1.22) \u(t,x)\ ^j'E&WixfV+'tpiiV + VVWv. 

Remark 1.4. As a byproduct of our method of proof, we obtain the limiting ab- 
sorption principle for the perturbed Dirac operator under assumption (1.21) (see 
Section 3.2). The LAP had been proved earlier for the free Dirac equation by Ya- 
mada [42] , and for the Dirac equation with potential (and with mass) in [28] under 
quite stronger assumptions. 
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2. The self-adjointness of the perturbed operators 

In this section we check the selfadjointness of the perturbed operators A\y and 
T>v under quite general assumptions on the potentials A, B, V, which in particular 
are implied by the assumptions of Theorems 1.1, 1.2 and 1.3. Most of the material 
here is standard; however we decided to include a sketch of the proof for the sake 
of completeness. Moreover, the use of Lorentz spaces techniques (see the Appendix 
for a short review) makes the proofs quite straightforward. 

It will be useful sometimes to express the magnetic laplacian both in the covariant 
form 

(2.1) H = -(V + iA(x)) 2 +B(x) 
and in the expanded form 

3 

(2.2) H = -A + W{x,D), W(x,D) =^2a j {x)d j +b(x) 

where 
(2.3) 

3 

aj(x) = -2iAj(x), b(x) = -i^djA^x) + \A(x)\ 2 + B(x), Aj,B e R. 
Then we have the following: 

Proposition 2.1. Consider the operator on Cg°(R") 

(2.4) H = -(V + iA{x)) 2 +B{x), 

where A(x) : R" — > R" and B{x) : R" — > R are measurable functions. Assume that 
the Lorentz (weak Lebesgue) norms of the coefficients 

(2.5) \\A\\ L n,oo < C , \\B\\ L n/2,oo<C 

are bounded by some constant Co > small enough. Then H has a (unique) self- 
adjoint extension to H 2 (R n ). 

Proof. Our proof is based on the standard results on quadratic forms, see e.g. the 
standard reference [29]. First of all we notice that by (2.5) we have immediately 

\A(x)\ 2 g £" Ao ° 

with a small norm. Now, the quadratic form q{4>, ip) given by 

q(<p, V>) = ((V + iA)<p, (V + iA)yj) L 2 + (Bp, ^P) L 2 

is well defined on the form domain H 1 under assumptions (2.5). Indeed, consider 
the identity 

(2.6) = IIVVIH2 + ((\A\ 2 +B)^) L 2 +23(AV^) L 2; 

using the embedding H 1 C x 2 "/(™- 2 ). 2 j the Holder inequality in Lorentz spaces (see 
the Appendix at the end of the paper for a quick synopsis of the relevant results), 
and recalling assumption (2.5), we have easily 

\qty,1>)\ <l|VV||i 2 + C\\\A\ 2 + B\\ Ln/ 2,oo\\Hi\\ L ^ + C\\A\\ Ln ,oo\\Vi> ■ 

<||v^||| 2 + cCoH\\ 2 L + cc 'o\\vM\l^M l ^,2 < c\m\\h. 

It is clear that the form is symmetric, since A and B are real valued. Now, recall- 
ing Theorem VIII. 15 in [29], in order to prove that q is the form associated to a 
(uniquely defined) self-adjoint operator, it will be sufficient to show that it is closed, 
i.e., its domain H 1 (M. n ) is complete under the norm 

(2-7) M 2 =QW>,1>) + C\Ml' 
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for some C > 0, and that it is semibounded, i.e., 
(2.8) qW,1>) > -C\\1>\\ 2 L i 

for some C > 0. Both properties follow from the identity (2.6); indeed, by estimat- 
ing as above we obtain easily 

q(M)> ||V^||£ 2 -CCo||VV||i 2 . 
In particular this implies that the norm (2.7) is equivalent to the H 1 (W n ) norm, 
provided Co is small enough, so that the form is closed; and this implies also that 

(2.8) is satisfied with C = 0. □ 

For the perturbed Dirac operator we have a similar result: 

Proposition 2.2. Let V{x) = V*(x) be a 4*4 complex valued matrix on R 3 . 
Assume that 

(2.9) |MUs,oo<C 0) 

for some Co > sufficiently small. Then the perturbed Dirac operator Vy = V + V 
is self-adjoint on if 1 (R 3 , C 4 ). 

Proof. The proof is analogous to the proof of Theorem 2.1. We define the quadratic 
form q : H 1 / 2 x H 1 / 2 — ► C associated to the operator Vy as 

V) := (p<p,il>) + (V<p,il>). 

First we prove that the domain of q is H 1 / 2 . With the same arguments of the 
previous theorem we estimate 

\q{<p,<p)\ < |M|^ 1/2 + cwvw^ooW^w^n^ 

< Mmn + c\\v\\ L ^y\\l 2nnn . 1) , 2 

< (i + c||y|U3,oc)||^|| H1/2 

(where we used the embedding H 1 / 2 C L 2 ™/(™ _1 ^ 2 ). From this point on, the proof 
proceeds exactly as in Proposition 2.1 □ 

3. The limiting absorption principle 

The essential tool in our proof will be the spectral theorem in the following ver- 
sion: given a sclfadjoint (unbounded) operator A on L 2 and a continuous bounded 
function /(A) on R, the operator f(A) can be defined as 

(3.1) f(A)4> = -- ■ L 2 - lim / f(X)QR(X + ie)<f>d\ 

IT ej.0 J 

for any G L 2 . Here R(z) = (A — z)^ 1 denotes the resolvent operator of A (see 
e.g. [38]). Under suitable assumptions on H, the limit operators R(X ± iO) = 
lim e |o i ie) are well defined as bounded operators in weighted L 2 spaces; this 
is usually called the limiting absorption principle (see below for details). Thus we 
have also the simpler representation 

(3.2) f(A)4> = ~ ■ J f(\)$3R(\ + iQ)(t>d\. 
Recalling the definition (6.3), consider now the operators 

3 

H = -A + W(x,D) = -A + Y,aj(x)dj+b(x) 

3 = 1 

and 

V v =V + V(x). 
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In Section 2 we proved that, under assumptions (2.5) on aj,b and V(x), both H 
and Vy are selfadjoint operators on L 2 . In particular, the spectral formula (3.1) 
holds for both. We shall use the following notations: the free resolvents will be 
written as 

R (z) = (-z- A)- 1 , R v {z) = {-zh + V)- 1 
while we shall use the notation R(z) for both perturbed resolvents: 
R(z) = (—z — A + W)-\ R(z) = (-z + V + V)- 1 . 

From the context the meaning of R(z) will always be clear. Note that Ro(z) is 
defined for all z £ K+ while Rv(z) is defined for z £ R, and the same properties 
hold for the perturbed resolvents. 

Our first task will be to show that the stronger representation (3.2), i.e., the 
limiting absorption principle, holds also for the perturbed operators. For A = —A 
this is a classical result (see e.g. Agmon [1]); here we shall use a very precise version 
of the principle, due to Barcelo, Ruiz and Vega [4]. On the other hand, for the Dirac 
operator only a few results are available, which concern the case with a nonzero 
mass term (see [28], [42]). 

The classical results on i? (see [1]) state that the limits 

(3.3) lim R (X±ie) = R (X±iO) 

exist in the norm of bounded operators from L 2 ((x) s dx) to H 2 ((x)~ s dx) for any 
s > 1; the convergence is uniform for A belonging to any compact subset of ]0, +oo[, 
and the following estimate holds 

(3.4) \\(x)- s R (\± l O)(x)- s f\\ L ,<^-\\f\\ L2 VAX), s>\. 

In n = 3 dimensions, the operators i?o(A ± iO) have the explicit representation 

1 f e ±iV\\x-y\ 

(3.5) R (X ± iO)g(x) = — / — r-g(y)dy, A > 0. 

47T J \x- y\ 

Recall also that for A < we have the similar formula 

1 f e^VT-^T \ x ~v\ 

(3.6) Ro(X)g(x) = —J _ g(y)dy, X < 0. 

These results were extended in [4] to more general weights. Introduce the norm 

(3.7) 1110(^)111 = sup / + - 2 k{r)r 2 , 1/2 dr where h(r) = sup \a(x)\. 

fJ.>0 J ft (T /i ) 1 \x\=r 

For any measurable function on R" such that supp/ C suppa, we can consider the 
(semi-)norm 

\\f\\ L 2 {a{x)dx) = \\a{xfl 2 f\\ L , < oo 

and we can define a Hilbert space L 2 (a(x)dx) as the closure in this norm of the sub- 
space of functions with support contained in suppa. Then we can summarize 
Theorems 1 and 2 in [4] as follows: 

Theorem 3.1 ([4]). Let a(x) be a nonnegative function on R™ with \\\a\\\ < oo, 
and denote by Rq(X ± iO) the limit operators (3.3). Then the operators Rq(z) for 
z ^ R + and R (X±iO) can be extended to bounded operators from L 2 (a(x)~ 1 dx) to 
L 2 (a(x)dx), and the following estimates hold: 

C 

(3.8) ||^o(A±iO)/|| L 2 (o(x)dx) < — =||H|| • \\f\\ L z(a(x)-idx), X^O 
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(here of course i?o(A ± iO) = i?o(A) for A < 0) 

(3.9) ||Vi2o(A±iO)/|| L2(o(x)tfa) < C|||o||| • \\f\\ L i{a{ x )-id x y 

Moreover, the limiting absorption principle holds in the weak form: for all f,g£ 
L 2 (a(x)- 1 dx) 

(3.10) lim(i? (A ± ie)f, g) = (R (X ± iO)f, g). 

Remark 3.1. It is not difficult to extend the estimates (3.8)and (3.9) to the whole 
complex plane. Indeed, fix two functions f,g e Cfi° with support contained in 
supp a and consider on the half plane 

S = {z: 3?z > 0} 

the holomorphic function 

(3.11) F{z) = z 1 ' 2 {R a {z)f,g). 

It is clear that F(z) is continuous on S up to the boundary, moreover it satisfies 
the estimate 

(3-12) |F(x)| < C||H|| • \\f\\ L i(a(x)-idx)\\g\\Li(a(x)-^dx) 

on the boundary 3z = 0, and finally it has a polynomial growth for \z\ — > +oo, 
as it easily follows from the explicit expression of Ro{z) as a convolution operator 
(see [4]). By the Phragmen-Lindelof Theorem (see e.g. [36]) on the half plane we 
immediately obtain that estimate (3.12) holds on all of S. A similar argument can 
be applied in the lower half plane 3z < 0. In conclusion we obtain 

C 

( 3 - 13 ) \\Ro(z)f\\^(a(x)dx) < —7= IIMH ' \\.f\\L 2 (a(x)- 1 dx) 



for all / € L 2 (a(x) 1 dx) (see also part (ii) in Theorem 1, [4]). Notice that this 
estimate holds on the whole complex plane, in the sense that we apply it to i?o(A± 
iO) when z e K+ . 

If we apply the same argument to the function 

G(z) = (VR (z)f,g) 

we obtain in an analogous way the estimate 

(3.14) ||Vi2o(«)/|| L 2 (o(x)tfa) <C|||o|||-||/|| ia(o(l) -i <fa) , zeC. 

We now specialize the theorem to a particular choice of weights. Precisely, 
consider the family of functions 

(3.15) wp(x) = |x|(|log|a;|| + l) /3 , [3 > 1. 
As it is proved in [4] (see Proposition 1), the norms 

II l 10 ^ 1 II I < +°° 

are finite for all [3 > 1, hence we can apply 3.1 with the choice 

a(x) = (wb(x))^ 1 = . — 1—r. —77. 

NI(|log|x|| + l)' 3 

In this case it is possible to improve the above result and to obtain a stronger 
version of the limiting absorption principle. To this end, we need the following 
Lemma, which is inspired by [1]: 

Lemma 3.2. Let H be a Hilbert space, H' its dual, and Hq a second Hilbert 
space compactly embedded in H' . Let Tj,T (j — 1,2, . . .) be bounded operators in 
C{H, H') such that 
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(i) Tj,T are symmetric for the pairing (•, -)h'xH, i-C-, 

{Tf,g) H 'xH = (Tg,f) H 'xH V/,g e H; 

(ii) Tj,T e C(H,Hq) and, for some constant C independent of j, 

\\Tj\\c { H,H ) < C 

Assume that 

(3.16) Tjf Tf weakly in H' for all f €H. 
Then Tj — > T in the operator norm of £(H,H'). 

Proof. Fix an / G H; the sequence Tj/ converges weakly to Tf in if', and is 
bounded in H by (ii), hence it admits a subsequence which converges in the norm 
of if', and the limit must be the same i.e. Tf. By applying the same argument to 
any subsequence of Tjf, we conclude that the entire sequence Tjf converges to Tf 
in the norm of if . 

Now, let fj be any sequence which converges to / weakly in H. Then we have 
for all g e H 

(Tjfj,g) = {Tjg,fj)^{Tf,g) 
since Tjg — > Tg strongly in if' and fj — / weakly in H . In other words, for any 
fj — ^ / weakly in if we have that Tjfj T/ weakly in if'. But, as in the first step, 
we can remark that the sequence Tjfj is bounded in Ho and by compact embedding 
we obtain that the convergence is strong: Tjfj — > T/ in the norm of i?'. 

By the same argument we obtain that, for any fj — ^ / weakly in if, the sequence 
T fj converges to Tf in the norm of H' . 

Finally, assume by contradiction that Tj does not converge to T in the operator 
norm of C(H,H'). This means that we can find a sequence fj € H with norm 
|| /j || if = 1 such that 

\\T j f j -Tf j \\ H ,>e>Q 
for some e independent of j. By extracting a subsequence we can assume that fj — i 
/ weakly in H, and by the above steps we immediately obtain a contradiction. □ 

Then we can prove: 

Proposition 3.3. Let wp{x), x e R" one of the radial weights (3.15) for some 
fixed (3 > 1. Then, for all A 7^ 0, the limits 

(3.17) \im Ro(X ± ie) = R„(\ ± iO) 

ej.0 

exist in the norm of bounded operators from L 2 (wp(x)dx) to H 2 (wp(x)~ 1 dx) and 
satisfy the estimates 

(3.18) p (A ± iO)f\\ L2{w -i dx) < \\fhi( W0dx ), VA ± 0, 

(3.19) l|Vflo(A±iO)/|| i2(V<jx) <C(b) \\f\\ L 2 {Wfjdx) . 

Proof. We apply Lemma 3.2 with the choices: H = L 2 '(w p(x)dx) , and hence H' = 
L 2 (wp(x)~ 1 dx) with the standard L 2 pairing; H — H 1 (wp (x)~ 1 dx) for some 
arbitrary (3 with (3 > (3o > 1; the norm of i?o of course is 

II/IIho = II^o 1/2 /II! 2 + II^ 1/2 v/II| 2 - 

Finally, as operators Tj we shall take (any subsequence of) the resolvent operators 
Ro(\ ± ie) as e I 0, while T = R Q (\± iO), for some fixed Ael. 

We now check the assumptions of the lemma. The compact embedding of H 
into H' is clear. Also the symmetry of the operators in the sense of (i) is evident. 
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The uniform bounds on Tj , T as bounded operators from H to H' are simply the 
estimates (3.13), (3.14) applied with the choice a(x) = wp(x)^ 1 . But it is clear 
that the estimate (3.13) implies also the following estimate 

(3.20) \\Ro(z)f\\ L 2 {w -^ x) < ^4 ll/I^W*)' yz * °> 

which is only apparently stronger, in view of the trivial embedding 

L 2 {wpdx) C L 2 (wf3 dx). 

In a similar way we have 

(3.21) \\VMz)f\\ LHw -i dx) < C(p ) \\f\\L*( Wpdx ). 

These inequalities show that assumption (ii) of the Lemma is satisfied. Finally, 
assumption (3.16) is nothing but the weak limiting absorption principle of Barcclo, 
Ruiz, Vega (see (3.10)). 

In conclusion, Lemma 3.2 implies that the limit (3.17) exists in the norm of 
bounded operators from L 2 (wpdx) to L 2 (w^ 1 dx) . Moreover, by the identity 

ARo(z) = -I- zRo(z) 

we obtain that the limit exists also in the norm of bounded operators from L 2 (wpdx) 
to H 2 (wp dx). The estimates (3.18) and (3.19) follow from the corresponding 
estimates for general z. □ 

3.1. The limiting absorption principle for the magnetic laplacian. In what 
follows, we shall focus on the case n = 3 exclusively. We follow the standard 
approach, based on the resolvent identity 

R(z) = (-z - A + W(x, D))- 1 = R (z)(I + WR (z))-\ 

Thus the main step of the proof will consist in inverting the operator / + WRo in 
suitable weighted spaces. We shall assume that the coefficients aj(x) and b(x) in 
W(x,D), defined as in (2.3), satisfy the assumptions 

(3 - 22) Mx)l - n(^(|iSni + i)^ - \x\*(\iog\x\\ + iy 

for some s £ [0, 1], (3 > 1 and some constant Co small enough. 
Our result is the following: 

Proposition 3.4. Assume the coefficients ofW(x,D) — J2 a j( x )^j + b(x) satisfy 
(2.3) (3.22) for some Co small enough, some s e [0, 1] and some (3 > 1. 

Then the operator I +W Ro is invertible on the weighted space L 2 {wp{x){x) 2s dx), 
and the inverse operators (I + WRo(z))^ 1 are uniformly bounded for all z G C. 
Moreover, the strong limiting absorption principle holds for R(z), in the following 
sense: 

(i) the boundary values 

(3.23) ]imR(X±ie)=R(X±i0) 

exist in the norm of bounded operators from L 2 (wp(x)dx) to H 2 (w^ s 1 (x)dx); 

(ii) the following estimate 

C(/3) 

( 3 - 24 ) \\R(z).f\\mw ls (x)dx) < —f= ■ \\f\\L*{w {x)-idx) 

holds for all z £ C, z^0. 
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Remark 3.2. In the case s = we recover exactly the strong limiting absorption 
principle proved in Proposition 3.3 above for the free operator i? - The additional 
weight (x) s was considered in view of the estimates that will be needed in the 
following section. 

Proof. Consider the operator 

W(x,D)R (z)f = ^a j (x)d j R (z)f + b(x)R (z)f; 

we estimate the two terms separately. 
First of all we have 

\\w 1 / 2 (x) s a j (x)d j R f\\ L ^ < \\w(3(x} s a 3 \\ L ^\\w^ 1/2 d j Rof\\ L 2 < CoWw 1 / 2 f\\ L 2 

by estimate (3.21), and this implies trivially 

(3.25) \\wl /2 (x) s a 3 (x)d 3 R f\\ L 2 < CoWw}/ 2 (x) s f\\ L 2. 

In order to estimate the electric term, we recall that, from the explicit expression 
of the free resolvent, we can write 



FT 1 ' 1 



Then we have 
(3.26) 

\\wl /2 b(x)R Q (z)f\\ L 2 < ||^ /2 fe(x)|| i2 ||i?o(z)/|| L =o < \\wl /2 b(x)\\ L 2-C 



Recalling Young and Holder inequalities in Lorentz spaces (see Theorems A. 2, A. 3), 
we have 



o*l/l 



< CH/IU^ = C\\w- 1/2 wl /2 f < C\\w-^\\ L e4w^f\\ L 2. 



1/2 | 



1/2 . 



Since w^ 2 e L 6 ' 2 for any (3 > 1 (Proposition A.4), (3.26) gives 

\\w}/ 2 b(x)R (z)f\\ L 2 < CWw^b^Ui ■ \\wl /2 f\\ L 2. 
Now, by assumption (3.22) on b(x) we have easily 

\\wl /2 b(x)\\ L 2 < CCo 

and we conclude that 

(3.27) ||^ /2 6(x)i? (z)/|| L2 < CCo ■ K'VlU*- 

In a similar way we have 
(3.28) 



w 1 J 2 {x)bR {z)f\\ L 2 < \\w^(x)b\\ L e\\Ro(z)f\\ L3 < \\w^{x)b\\ L& -C 



1/2 



1/2, 



l/l 



L 3 



and 



R * 1/1 



L 3 



<C\\f\\ L * =C\\w- 1/2 (x)-'w}/ 2 (x)f\\ L1 



As above, we notice that w 



-1/2/ 



1 e L? for any (3 > 1, hence we have from (3.28) 



\\wl /2 (x)bRo(z)f\\ L 2 < C\\w 1 J 2 {x)b\\ L6 ■ \\w}/ 2 (x)f\\ L 2. 
Assumption (3.22) guarantees that 

\\wl /2 (x)b(x)\\ L e<CC 
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and, in conclusion, 

(3.29) \\(x)w}/ 2 b(x)Ro(z)f\\ L 2 < CC ■ \\(x)w}/ 2 f\\ L 2 

If we interpolate between (3.27) and (3.29), we obtain the estimate 

(3.30) \\(x) s w}/ 2 b(x)R (z)f\\ L 2 < CC ■ \\(x) s w}/ 2 f\\ L 2 
Summing up, from estimates (3.25) and (3.30) we get for all z € C 

(3.31) \\(x) s w l /2 WRo(z)f\\ L 2 < CC ■ ||(z) s ^ /2 /IU- 

Then it is clear that we can invert the operator / + WRo by a Neumann series on 
the space L 2 ((x) 2s wpdx). Hence, the standard representation 

(3.32) R(z) = R (z)(I + WR (z))- 1 

is valid. To conclude the proof of the Proposition, it is now sufficient to remark 
that, from property (3.17) of Proposition 3.3 and the uniform bounds on the norm 
of (/ + WRn(z))^ 1 we have just obtained (for s = 0), the limits in 3.23 exist in a 
weak sense. Proceeding as in the proof of Proposition 3.3, using Lemma 3.2, we 
deduce (i). Finally, (ii) is a consequence of (3.32) and the corresponding estimate 
(3.20) for R . □ 

Remark 3.3. Note that the assumptions of the preceding proposition can be ex- 
pressed in terms of the original coefficients A, B as follows: 

(3-33) \A{x)\ < — — -g^ ■ — , \VA(x)\ + \B(x)\ < ( '" 



a;|<a;>*(|log|z|| + l)0' 1 Wl 1 Wl " \x\ 2 (\ log |x|| + 1)0 
for some [3 > 1 and a constant Co > small enough. 

3.2. The limiting absorption principle for the Dirac operator and its per- 
turbation. In this section we will study the limiting absorption principle for the 
massless Dirac operator V; this property was studied by Yamada in [42] for the 
operator with mass. Moreover, as in the case of the magnetic Laplacian, we will 
extend this result to the perturbed operator T>y = V + V(x), under a suitable 
assumption on the potential V. 

It is well known that the spectrum of the free operator V is the whole real line. 
Due to the relation V 2 = —AI 4 , we immediately obtain the representation 

(3.34) R v (z) = R (z 2 )(V + zh), 

for all z G C with ?Rz = 0. Using this formula and the Proposition 3.3, we easily 
prove the following: 

Proposition 3.5. Let wp{x), i£R 3 be defined as in (3.15), for some fixed (3 > 1. 
Then, for all AeR, the limits 

(3.35) limi?x>(A±ie) = R v (\±i0) := R (X 2 ±i0){V + XI 4 ) 

ej.0 

exist in the norm of bounded operators from L 2 {wp{x)dx) to H 1 (wp(x)~ 1 dx) and 
satisfy the estimate 

(3-36) \\Rv(z)f\\ L 2^ w , j(x yi dx ^ < \\f\\L*(wfi(x)dx)> 

for all z e C. Moreover, we have the explicit representation 

R„(\±io)f = ^ [ 4r^? \h - »iM \m d v 




(3-37) ^— - ^>> J ^f/( 2 /)dy.. 
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Proof. The strong convergence of Rj>(\±ie) to Rj>(\±iO) in the space of bounded 
operators from L 2 (wf3{x)dx) to H 1 (w/3(x)~ 1 dx) is obtained by interpolation us- 
ing the property (3.17) and the representation (3.34); estimate (3.36) immediately 
follows from (3.34) and the estimates (3.18), (3.19), (3.20), (3.21). In conclusion, 
recalling the explicit representation (3.5) for Rq(X ± iO), after an integration by 
parts we get the formula (3.37) and this concludes the proof. □ 

At this point, we will proceed in a similar way to the case of the perturbed 
Laplacian and we will prove that it is possible to extend the above result to small 
electric perturbations of the free Dirac operator. As for the magnetic coefficients 
of W(x, D), we need to assume that the potential V satisfies 

(3-38) \V(x)\ < I i / \ n , C °> n — tttti 

|a;|(a;) s (|log|a;|| + l)* 3 ' 

for some s £ [0, 1], [3 > 1 and some constant Co small enough. We prove the 
following result: 

Proposition 3.6. Assume the potential V satisfies (3.38) for some Co sufficiently 
small, some s e [0, 1] and some [3 > 1. 

Then the operator I +VRj> is invertible on the weighted space L 2 (wp(x)(x) 2a dx) , 
and the inverse operators (I + VRv(z))^ 1 are uniformly bounded for all z e C. 
Moreover, the strong limiting absorption principle holds for R(z), in the following 
sense: 

(i) the limits 

(3.39) lim R(X±ie) = R(X±iO) 

exist in the norm of bounded operators from L 2 (wf3(x)dx) to H 1 (w^ 1 (x)dx); 

(ii) the following estimate 

(3-40) f\\L 2 (w (x)- l dx) < ■ \\f\\L 2 (w (x)dx) 

holds for all z € C, z^O. 

Proof. The argument is the same of the proof of Proposition 3.4 for the magnetic 
part of W. First we observe that, by hypothesis (3.38), we have 

\\w}/ 2 (x) s V(x)R v f\\ L2 < |K(x} s F(x)|| L HI V /2i? ^lli 2 ^ C ° ' \H 1/2 fhi- 

Hence we obtain the estimate 

\\w 1 J 2 {xYV(x)R T) {z)f\\ L , < \\w}J 2 (xyf\\ L2 , 

uniformly in z e C; thus we can invert the operator I + VRd by a Neumann series 
on the space L 2 (wpdx). Again, we can exploit the representation 

(3.41) R(z) = R v (z)(I + VR v (z))- 1 . 

By property (3.35) of Proposition 3.5 and the uniform bounds of (I + VR-p)^ 1 , it 
follows that the limits in (3.39) exist in a weak sense. Then we can procede as in 
the previous cases, using Lemma 3.2 and obtain (i). In conclusion, the estimate (ii) 
is an immediate consequence of (3.41) and the inequality (3.36). This concludes 
the proof. □ 

In the following we shall also need a weaker version of the last result: we shall 
require that V satisfies 

(3 - 42) " ^^(^(llog^ll +l)/3/2' 

for some s > h, f3 > 1 and some constant Co small enough. Then we have 
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Corollary 3.7. Assume the potential V satisfies (3.22) for some Co sufficiently 
small, s > \ and [3 > 1 . 

Then the operators I + VR-p are invertible on the space L 2 {{x) 2s alx) , and the 
inverse operators (I + V R-uiz))^ 1 are uniformly bounded for all z e C. Moreover, 
the strong limiting absorption principle holds for R(z), in the following sense: 

(i) the limits 

(3.43) lim R(X±ie) = R(X±iO) 

exist in the norm of bounded operators from L 2 ((x) 2s dx) to H 1 ((x)~ 2s dx); 

(ii) the following estimate 

(3-44) \\R{z)f\\ L i{(x)-^dx) < C ■ ||/||i,2(( x )2 S(te) 

holds for all z € C, z ^ 0. 

Proof. The proof is analogous to the proof of Proposition 3.6. Indeed, from estimate 
(3.36) and assumption (3.42) we have immediately 

\\(xyvR v \\ L2 < \\(xyw}/ 2 v\\ L ~\\w- 1/2 R v f\\ L2 < C \\wl /2 f\\ L2 

and by the trivial inequality 

w' 2 < C s (x) s , 

valid for all s > 1/2, we conclude that 

\\{x) s VRv\\ L 2<C \\{x) s f\\ L2 . 

Thus we can again invert (I+VR-p) with a Neumann series, and proceeding exactly 
as before we obtain the proof of the Corollary. □ 

4. Resolvent Estimates 

In this section we prepare the crucial resolvent estimates that will be used in the 
proof of the main results. In order to use the spectral formula, we need estimates on 
the perturbed resolvent operators and their derivatives with respect to A as bounded 
operators from suitable weighted LP spaces to L°° . We shall use the Holder and 
Young inequalities in Lorentz spaces extensively; for the convenience of the reader, 
we give a sketch of the main usefule results in the Appendix A. 

We consider first the resolvent of the magnetic laplacian. We recall that, by 
Proposition 3.4, the operators R(X±iO) = R (X±iO)(I+ W(x, D)R (X±iO)y 1 are 
well defined as bounded operators from L 2 (wp(x)dx) to H 2 (w p{x) _1 dx); moreover, 
we have the explicit representation (3.5). Our first result is the following: 

Lemma 4.1. Let R(X±iO) = R (X±iO)(L+W(x, D)R (X±iO))- 1 be the resolvent of 
— A + W and assume the coefficients ofW(x,D) = a j( x )dj + b(x) satisfy (3.22). 
Then, for all X > 0, the following estimates hold: 

(4.1) H^AiiOj/llio. <C\\w}/ 2 f\\ L 2, 

(4.2) ||0 A i2(A±iO)/|| L ~ <c(l + -^=) \\(x)w}/ 2 f\\ L2 . 

Proof. The estimate (4.1) is the easiest one. In fact, by formula (3.32) and the 
explicit representation (3.5) for i? , we obtain 

||i?(A±*0)/|| L oo < C- H^* |(7 + WRo)- 1 /^; 
\x\ 
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using Young inequality in Lorentz spaces, we get 

iii?(A±io)/i| L oo < wd+wRor'fhs^ 

The uniform bound for the operators (7 + WRo)^ 1 proved in Proposition 3.4 and 
the observation that wj^ 2 G L 6 ' 2 , for all (3 > 1 (see Proposition A. 4) are sufficient 
now to conclude the proof of estimate (4.1). 

In order to proceed with the proof of (4.2) we observe that from (3.5) we imme- 
diately obtain the following explicit representations, for all A > 0: 

(4.3) d x R (X ± i0)f = R 2 (X ± i0)f = ±^= J™ e ± ^ x -y\f(y)dy, 

(4.4) d 3 Rl{\ ± i0)f = ±^l°° e*^*-* J2 ^lA f{y)dy - 
At this point, differentiating in (3.32) we get 

(4.5) d x R{\ ± iO) = A + B 
where 

A = R 2 (X ± iO)(I + WR (X ± iO))- 1 

and 

B = R (X ± *0)(J + WRq(X ± *0)) _1 Wi2§(A ± *0)(J + VUR (A ± *0)) _1 . 
We treat separately the two terms. By (4.3), we estimate 

p/ikoo < ^i+wR^nw 

< -^||( a ;)- 1 U ; /3 ( a; )- 1 /2|| i2 ||( x ) U ; /3 ( x )i/2 (/ + ^ o) -i / || L2 . 

We observe (Proposition A. 4) that (x)~ 1 W/}(x)~ 1 ^ 2 G L 2 for all (3 > 1 and, by the 
uniform bound for the norms of (I + WRo)~ 1 in the space of bounded operators 
onto L 2 {{x)wp{x)dx) for (see Proposition 3.4), we conclude that, for some C > 

(4-6) ||A/|| L oo < ^=\\{x)w ( X y/ 2 f\\ L2 . 

For the estimate of the term B, we start with some computation on the operator 
WR 2 . Using the representation (4.4), we obtain 

Ww^a^RlfU* < \\wy\i\\LA\diI%f\\L~ < C- Ik^lUHI/IUi. 
By the above observation that 

||/|| Ll <||(x)^ 2 (x)/|| L2 , 
it turns out that, if w X J 2 aj G L 2 , then 

(4.7) \\w fj (xy/ 2 a,(x)d,R 2 f\\ L2 <C-\\(x)wl /2 (x)f\\ L2 . 

In a similar way, using (4.3), we have 

1/2^2,11 . / II l/2w mi „ D 2 ... ^ C 1/2, 



\\w^bR 2 f\\ L , < \\w^b(x)\\ L 4R 2 f\\ L oo < — \\ w ^b\\ L2 \\f\\ L1 . 
If we assume that w\ b £ L 2 , we conclude that 
(4.8) \\w (xy/ 2 b(x)R 2 o f\\ L2 < ■ \\(x) W f(x)f\\ L2 . 
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Inequalities (4.7) and (4.8) can be unified now, to show that, under the assumptions 

(4.9) w^ajeL 2 , wl /2 beL 2 , 
the estimate 

(4.10) \\w (x)^ 2 W(x,D)R 2 (X±iO)f\\ L ,<c(l + ^j \\{x)wp{xfl 2 f\\ L , 

holds, for some C > 0. Observe that assumptions (4.9) are weaker than (3.22), so 
that they are obviously satisfied by the hypothesis of the Lemma. 

Now we are ready for the estimate of the term B. First, we use the representation 
(3.5) for Ro to obtain 

< ^-^[(I + WRo^WRKl + WRo)- 1 ^ 

< \\{I + WRv)- 1 WRl(I + WRV)- 1 =: \\Tf\\ L3/2A . 
As before, we use the properties of the weights wp(x) to observe that 

||.g|| L 3/2,i < \\w fj {x) 1/2 g\\ L 2. 
Then, the last series of inequalities gives 

< \\w (x)^ 2 Tf\\ L 2. 

Now we use the uniform bounds for the inverse operators (/ + WRq)^ 1 (see Propo- 
sition 3.4) to proceed with 

< Ww^x^WRKl + WR )- 1 f\\ L 2-, 

finally, by inequality (4.10) and the above mentioned estimates on the norms of 
(1 + WRq)^ 1 in the space of bounded operators onto L 2 ((x)wfj(x) 1 ^ 2 dx), we obtain 
the estimate 

(4.11) ||B/|| L « <c(l + -±=^ \\{x)wp{xf' 2 f\\ L ,. 

In conclusion, estimates (4.6), (4.11) and the representation (4.5) conclude the proof 
of (4.2) and the Lemma. □ 

Remark 4.1. The limiting absorption principle allows us to rewrite the spectral 
formula in the following way: for any (smooth, compactly supported) function 
<t>{\) on K, and any test function /, 

r+oo 

(4.12) (f>(-A + W)f= / <t>{X)%R(\ + iO)fd\. 

Jo 

where the integral is restricted to the positive real axis since of course QR(X) = 
for negative A. 

The resolvent estimates just proved imply that we can integrate by parts in the 
above formula, i.e., if 

</»(A)=V'(A) 

then 

r+oo 

(4.13) <j>(-A + W)f = V'(A)3?i?(A + i0)/dA 

Jo 

ip{X)d x ^R{X + i0)fd\ 



f 

Jo 
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The problems arising from the singularity at A = are easily overcome. To prove 
this, consider a cutoff function x(A) supported in [—L,L], and write 

,+oo 

<f>(-A + W)f= lim / (f>(X)(l-x(\L))QR(X + iO)fd\ 

whence 

,1/L 

4>(-A + W).f = - lim L V(A)x'(Ai)3?i?(A + iO)/dA 

L^+oo J 
,+oo 

- lim / (l-x(\L))iP(\)d x QR(\ + iO)fd\. 

L^+oo J Q 

=u L + v L . 

The last term vl converges to (4.13) uniformly, thanks to estimate (4.2) (and 
Lebesgue's dominated convergence theorem), hence it is clear that tit = <f>(— A + 
W)f — vl also converges uniformly, and it will be sufficient to show that its limit 
is 0, e.g., in distribution sense. To estimate the integral 

,1/L 

u L = -L / V(A)x / (AL)Q J R(A + iO)fdX 
Jo 

we can use the identity 

(4.14) 3?i?(A + i0) = (/ + i?o(A - i0)VK)" 1 3?i? (A + iO)(J + WR {\ + i0)) _1 . 
Consider then the L 2 product 

(9fi2(A + iO)/, <?) - (3i? (A + »0)(J + VFi? (A + iO))" 1 /, (/ + ^o(A + tO))- 1 *?). 
From the explicit formula 

Si? (A + ^ = c/^^ %)rfy 
J \x-y\ 

we have 

|3i?o(A + i0)h| <C\/A y |%)|cfy 

which implies 

||3?i? (A + i0)fr|| L oo < C\/A||/i|| L i <C>/A||<a;)tuJ /2 /i|| L 2 

for any f3 > 1. Recalling now the uniform bound for (/ + W^Ro(A + «0)) _1 in 
Proposition 3.4 in the weighted L 2 norms with weight (x)w 1 J 2 , we obtain easily 

\(<ZR(X + iO)f,g)\ < CVXMx^fhqWw^gh,. 
From this estimate it is easy to prove that 

,1/L 

(u L , g) = -L ^(A) x '(AL)(9i?(A + iO)/, fl )dA -> 
Jo 

as i — > +oo, which concludes the argument. 

We will prove now an analogue of Lemma 4.1 for the Dirac operator. In what 
follows, R(z) = {—zli + T> + V)^ 1 denotes the resolvent of the perturbed Dirac 
operator. Our approach here will be slightly different: we shall use the formula 

(4.15) R(z) = R v (z) + R v {z)V(x)R v (z){I + V(x)R v {z))-\ 

valid for all z € C (to be interpreted of course, for z = A e R, as the extended 
resolvents R(X) := R(X ± iO) on the weighted L 2 spaces, as given by Proposition 
3.6 and Corollary 3.7). When inserted in the spectral formula, the first term Rx> 
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at the right hand side reproduces the solution to the free Dirac equation, and the 
main part of our proof will be the estimate of second term 

(4.16) Q~R v VR v {I + VR v )-\ 

To this end, we shall need an explicit representation for Rj>(X±iO), which is easily 
obtained from the formula 

(4.17) Rt>(\ ± iO) = i? (A 2 ± iO)(V + Xh). 
Recalling (3.5), after an integration by parts we obtain 



iX f e ±lX \ x -y\ I T ^ - y 3 



4tt J M 3 \x-y\ I jri \x-y\ ' 

i f p ±i\\x-y\ 3 — 

(4-18) +i- / is-E^-r— T^) d f- 

4?r J K 3 |a;-y| 2 ^— J \x-y\ 

From here we derive immediately an analogous representation for 

d 

Rt>W = -q^Rv(X); 
indeed, differentiating (4.18) with respect to A, we get 

R%{X±iQ)f = Lj 4 + £ a . 




f(y)dy 



(4.19) 

We collect all the necessary estimates in the following lemma (we write for simplicity 
i?x> (A) instead of Rd (A ± iO) since the estimates are the same) : 

Lemma 4.2. Suppose that 

Cn 



(4.20) \V(x)\ < 



x^ixyQiogixW + iy/ 2 ' 



for some s > |, /?> 1, Co > 0. Then the following estimates hold for all e > 
small enough and all A € R: 

(4.21) || (a:) 1//2+£ yi?f,(A)/|| L 2 < C £ • (A) • || (*} 3/2+ 7IIl 2 , 

(4.22) ||i?x>(A)V r i?x'(A)/||L=c < C £ • (A) 2 • \\(x)^ 2+ \f\\ L2 

(4.23) ||i&(A)Wfc(A)/|| L ~ + ||ifc(A)W&(A)/|| i0 . < C t • (A) 2 • || {x)^f\\ L , 
for some C = C e independent of A. 

Proof. In the following we shall use the shorthand notation, for s € R, 
(4-24) \\f\\L>:=\\(x)V\\ L ' 

From the explicit representations (4.18) and (4.19) we have the simple pointwise 
estimates 

(4.25) \Rv(X)f\ < C(|A| • M" 1 + |z|- 2 ) * /, \R^(X)f\ < C(\X\ + |a:| _1 ) * /. 
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Since \x\ 1 <G L 3 -°° , by the Young inequality in Lorentz spaces (see the Appendix) 
we get 

wvRiwfh* < \\v\\ L , ■ iai ■ hi * /ik- + immikr 1 * /ik- 

<l|V|k»(|A|- H/Iki + 11/11^,0. 
By the obvious inequalities valid for all e > 

(4.26) ||/|| L1 < C(e)||/||^ /2+e , H/lkv^x < C(e)||/|| iW , 
we arrive at the first estimate 

(4.27) l|ra|>(A)/|k, < C(e)||y|k,(A)||/|| Li/2+e . 

Since || V|k 2 < oo by assumption (4.20) as soon as 7 = 1/2 + e < s — 1, we see that 
(4.21) follows provided e is suitably small. 

In a similar way, in order to prove (4.22) we use again (4.25) and we write (recall 
that |.t|- 2 g L 3 / 2 <°°) 

\\Rv(\)VR v (\).f\\L~ <C(\X\ ■ \Wx\~ 1 *VRvf\\ L oo + \\\x\~ 2 *VRvf\\ L oo) 

<C(\X\ ■ \\VR v (X)f\\ L s,^ + ||Wfc(A)/|ks,i) . 

For the first term we can write, recalling again (4.25), 

(4.28) \\VR v (X)f\\ L3/ ^ < ||y|| i3/2 ,i|A| • ll^r 1 * f\\ L oo + \\V\\ L 2\\\x\- 2 * f\\ L e, 3 

< ll^lks/^IAI-ll/lks/^ + ll^ll^ll/ll^ 

^(IIVlks/^IAl + HVlkOII/lk^. 
(see (4.26)), while for the second term we have 

(4.29) ||Wfc(A)/lk3.i < ||V|k..i|A| • UN" 1 * /Ikoc + \\V\\ L s, 2 \\\x\-2 * f\\ L8 , a 

<ll^|k..i|A|-||/|ks/ a .i + ||V|ka. 2 ||/|k a 
<(||V|ks. 1 |A| + ||V|ke, a )||/|k S/2+< 

where we have used (4.26) and the trivial inequality ||/|k 2 < ll/IU?,, V7 > 0. 
Summing up, we get 

(4.30) \\Rv(\)VR v (\)f\\ L ~ < C ■ C(T0<A) 2 ||/|| Li/2+e 
where the quantity 

(4.31) C(V) := Hylks/2,! + ||V|ks,i + ||V|k«, a + ||y|| L 2 < 00 

is finite by assumption (4.20) (see also the Appendix A). 
The proof of (4.23) is similar: by (4.25) we get 

\\Rl(\)VR v (\).f\\L~ < C (|A| -||1 * VRdSWl- + \Wx\- 1 * VR v .f\\ L ™) 

<C(\X\ ■ \\VR v (X)f\\ Ll + ||Wfo(A)/|ks/ 2 ,i) . 

We have already estimated the second term in (4.28), and for the first one we have 

(4.32) \\VRv(\)f\\ L i < \\V\\ L ,„\\\ ■ \Wx\- 1 * f\\ L3 + mksllkr 2 * /Iks/2 

<(||^lk3/ 2 |A| + ||V|ks)||/|k. 
<(ll^lks/ 2 |A| + ||V|ks)||/|k w< 

and hence 

(4.33) \\Rl(\)VRv(\)f\\ L ~ < C ■ C'(y)(A) 2 ||/|| L2/2+£ 
where the quantity 

(4.34) C'(V) := \\V\\ La/ 2 + \\V\\ La/ 2.i + \\V\\ L s + \\V\\ L 2 < 00 
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is finite again by assumption (4.20). 

Finally, the last estimate can be obtained as follows: 

\\R v (X)VR 2 v (X)f\\ L ~ <C(\X\ ■ \Wx\- 1 * VRlf\\ L oo + \\\x\- 2 * VRlf\\ L oo) 

<C(\X\ ■ WVRvWfhw + \\VR v (\)f\\ L3 ,i). 

Proceeding as above, we estimate 

(4.35) ||W&(A)/|| L 3/ a ,i < IIVIUs/^iIAI -||1 * f\\ L ~ + liyll^^HH- 1 * f\\ L oo 

< 11^11^4(^(11/11^ + 11/11^,0 

< ||^|| i ./ 2 ,i4<A)||/|| iS/2+< 

and 

(4.36) \\VR 2 v {X)f\\ L ^ < ||V|| i ,, 1 |A|-||l*/|| i oo + ||V|| i 3,i||| a: |- 1 */|| L cc 

< 11^11^x4(^(11/11^ + 11/11^,0 

< ll^|Us.i4<A)||/|| LWe 

whence 

(4.37) ||ifc(A)W&(A)/|U~ < C ■ C"(y)(A) 2 ||/|| L , /2+£ 
where the quantity 

(4.38) C"(V):=\\V\\ L ,,2,i + \\V\\ L *.i < oo 

is finite by assumption (4.20). □ 

Remark 4.2. The same remark concerning the simpler version of the spectral for- 
mula (4.12) and the integration by parts formula (4.13) applies also to the Dirac 
resolvent, with obvious modifications in the proof. 

5. Proof of Theorem 1.1 

Let (<Pj)j_ 1 be a standard Paley-Littlewood partition of the unity, with the 
properties 

(5.1) Vj(A) = <£o(2~ J A), Po + J>i = l, 

for a suitable tpo £ C5 30 . We consider the Cauchy problem 

f u tt (t, x) - Au(t, x) + W(x, D )u = 

\ u(0,x) = 0, u t {0,x) = <pj(y/-A + W)g(x), 

The solution can be represented using the spectral formula as follows: 



1 f + °° 

(5.3) u{t,x) = —J Vj (VX) 



sin(iVA) 



R(X)gdX, 



VA 

and after an integration by parts (see Remark 4.1) this gives 

(5.4) u(t,x) = — cos(iVA) dx<Pj(V\)R(X)g + <Pj(V\)dxR(X)g dX. 

t Jo L - 1 

Thus, recalling estimates (4.1) and (4.2), we have 

\u(t,x)\ < ^\\( X } W l/ 2 g\\ L2 (\d X < Pi {V\)\+ + \<Pi{VX)\^dX 

and a change of variables A = 2 2j ^ in the integral gives immediately 

(5.5) \u(t,x)\<j2^\\(x)wfg\\ L2 
with some constant C independent of j and g. 
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If we now define as usual 

<Pj = <Pj-i + <Pj + Vj+i, <P-i = 0, 
so that ipj = tpjipj, we see that the Cauchy problem (5.2) can be written equivalently 



(5.6) 



J u tt (t,x) - Au(t,x) + W(x,D)u = 
I u(0,x) = 0, u t (0,x) = ip j {y/-A w )!p j (y/-A w )g(x), 

hence our estimate (5.5) implies also the estimate 

(5.7) \u(t,x)\ < j2^\\(x)w 1 f) /2 ^( x rA^)g\\ L2 . 

Finally, consider the original Cauchy problem (1.6), and decompose g as a sum 
g = J2j>o fjW~^w)g(x). By estimate (5.7) we obtain easily estimate (1.9). 

(5.8) \u(t,x)\ < ^2^||( 2 ;)^ /2 ^( v ^A^). g || L2 . 

j>0 

The computations in the case of initial data of the form 

u(0,aO = /, u t (0,x) = 

are completely analogous, and we thus obtain estimate (1.12). 

Remark 5.1. In view of the application to the Dirac system, the following remark 
will be useful. If the initial datum g has the form 

(5.9) g = (-A w ) s h 

for some s > 0, a direct application of estimate (5.8) would give only 

(5.10) |«(t,aO| < y^2^||(x)^ /2 ^(v /Z A^)(-A w ) s / l || i2 . 

Actually, if we go back to the spectral formula (5.4), we see that the solution can 
be written 

Q r+oo 

(5.11) u(t,x) = — A s/2 cos(n/A) dx<pj(V\)R(\)h + <pj(V\)d\R(\)h dX. 

t Jo L J 

with an additional factor X s / 2 . Thus, proceeding as above, we arrive at the simpler 
estimate 

(5.12) \u(t,x)\ Kj^^mw^^i^TA^M^. 

j>0 

We now prove estimate (1.11) under the stronger assumption (1.10) on the po- 
tential W(x,D). Consider first the case of initial data of the form 

u(0, x) = 0, ut(0, x) = g. 

We can write g as follows: 

g = (1 - A + W)- 1 -^ - A + W) 1+f g 
for some fixed e > 0. Then the solution u can be represented as 

u (t iX ) = ±- [ + °° ^{Vx f^ 1 ^ R(X)hdX 

where 

h = (1 - A + W) 1+e g, tP(VX) = (l + A) 1+e . 
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Proceeding as above, after an integration by parts we arrive at 

s~i r+oo 

\u(t, x )\<^\\( x )wl /2 h\\ L2 / ((l + A)- 1 - e + (l + A)- 2 - £ )dA 
£ Jo 

and hence 
(5.13) 

\u(t,x)\ < j\\(x)wl /2 (l-A + W) 1+ <g\\ L . < j\\(x)V 2+ *(l-A + W) 1+ <g\\ L2 . 
To conclude the proof of the Theorem, it remains to show that 
(5.14) \\(x) 3 ' 2+e (l-A + W) 1+e g\\ L , <\\{xfl 2+ *g\\ H ^. 

We start from the inequality 

\\(xni-A + w)f\\ L2 <\\( x yf\\ H2 

which is obviously valid for any s > 0. By a standard complex interpolation argu- 
ment, interpolating with the trivial inequality 

||<*>VIU* < ||<a>VIU* 

we obtain that 

IK^'Cl - A + W) e /IU=» < ||<o:>-/ll^. 
for all < e < 1 and all s > 0. This implies 
(5.15) 

\\(x)'(i-*+w) 1+e f\\ L , < \\(x) s (i-A+w)f\\ H ^ < \\(xyf\\ H2+2e +\\( x ywf\\ H 2,. 

The last term is of the form 

(5.16) \\(x) s W(x,D)f\\ H 2, < \\(x) s a(x)Df\\ H 2 C + ||<a:)'6(a:)/|| ff2 .; 
in order to estimate it, we recall the Kato-Ponce inequality (see [21]) 

(5.17) \\(D)«(vw)\\lp <C\\(Dyv\\ Ln \\w\\ LP2 + C\\v\\ L p 3 \\(D)i W \\ LP4 

which is valid for all q > 0, jT x — p\ + p^ 1 = p^ 1 + p^ 1 . With the choices 
v(x) = a(x), w{x) = (x) s Df(x), q = 2e, p\ = pz = oo and p 2 = P4 = 2, we obtain 

\\{D)*{x)'a(x)Df\\„ <C\\(D) 2 *a\\ L °o || {xYDf\\ L .+C\\a\\ L oo\\{D) 2 *{{ X yDf)\\ L .. 

Now it is clear that 

\\{D)^{{xyDf)\\ L ,<c\\{xyf\\ H ^ 

(use again complex interpolation between the cases e = and e = 1) and in con- 
clusion we obtain 

\\(D)*(x)'a{x)Df\\ L , < C\\(D) 2t a\\ L - \\{x)'f\\ m+ ^. 

Here we have used the simple fact that 

NU~ <q|(^> 2e a|U~. 

The corresponding estimate for the electric term is analogous (actually simpler): 

\\{D) 2e (x)'b{x)f\\ L , < C\\(D) 2 *b\\ L °o\\(xyf\\ H 2<. 

Recalling now (5.15) and (5.16) we conclude the proof of estimate (1.11). 
On the other hand, when the data arc of the form 

u(0,x) = f, u t (0,x) = 

the computations are completely analogous and we obtain estimate (1.14) under 
the stronger assumptions (1.13) on the coefficients. 
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6. Proof of Theorem 1.2 

Remark 6.1. We notice that Theorem 1.1 (and Remark 1.3) can be trivially ex- 
tended to a system of wave equations of the form 

(6.1) u tt - (V + iA(x)) 2 u + B(x)u = 

where u(t,x) is a C N valued function and A\(x), A 2 {x), A 3 (x), B{x) are C NxN 
matrices whose coefficients satisfy the assumptions of the Theorem. The resulting 
dispersive estimates have exactly the same form as in the scalar case. 



Consider now the Cauchy problem 
(6.2) 



iu t — Vu — V(x)u = 
^ u(0,aO = /(a;). 

If we apply to the pertubed Dirac system the operator id t + V + V we obtain that 
u is also a solution of a 4x4 system of perturbed wave equations of the form (6.1) 
with 

(6.3) A j (x) = ~(a j V(x) + V(x)a j ), 

(6.4) B(x) = VV{x) + V{xf + Aj + A 2 2 + A 2 3 + ®i A 3 
and initial data 

(6.5) u(0,x) = f, u t (0,x)=r 1 (V + V)f. 
Note that the perturbed operator 

(6.6) -A w = -(V + iA(x)f +B(x) 
is exactly the square of the operator V + V: 

(6.7) -A W = (V + V) 2 
and hence the initial data for (6.1) can be written 

(6.8) u(0,x) = f, u t (0,x)=i- 1 (-A w ) 1 ^f. 

We are in position to apply to the solution u the estimates already proved in 
Theorem 1.1; keeping Remark 5.1 into account, we arrive easily at the estimate 

(6.9) \u{t,x)\ < jJ2^\\(^l /2 ^ + V)f\\ L ^ 

j>0 

provided the coefficients aj(x) and b(x) satisfy the assumptions (1.8). Recalling the 
explicit form (6.3) of the coefficients in terms of V(x), we see that V must satisfy 
the conditions 

]V[ n ~ \x\(x){\\og\x\\ + \)P 
from the magnetic term, and 

\ 2 1 i I r.T^/„\l ^ Co 



\V{x) 2 \ + \DV{x)\ < 



I^PdloglxH + l)/^ 

from the electric term, for some (3 > 1 and some small constant Co- Summing up, 
we obtain that (6.9) holds under assumption (1.17). 

The estimate in terms of the Sobolev norm can be obtained in exactly the same 
way as for the perturbed wave equation. Indeed, proceeding as in (5.13) we arrive 
at the estimate 

(6.10) | u (i,aO| < j\\(xf/ 2 +*(-A w )V 2 +*f\\ L2 . 
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The same arguments used at the end of Section 5 give here 
(6.11) \u(t,x)\<j\\(x)V 2 +*f\\ H3+2 < 
provided 

{D) 1+2e Aj e L°°, (D) 1+2e B e L°°, 

which is implied by 

(D) 2+2e V € L°°. 
7. Proof of Theorem 1.3 

By exploiting the connection between the massless Dirac and the wave equation, 
it is easy to obtain an optimal dispersive estimate in the unperturbed case. Indeed, 
let u(t, x) be a smooth solution of the free massless Dirac equation 

(7.1) iut(t, x) = Vu(t, x) 
with initial data 

(7.2) u(0,x) = f(x). 
Recall now the identity 

(id t +V)(id t -V) = (A-9 t 2 t )/ 4 ; 

if we apply the operator idt + T> to the system (7.1) we immediately obtain that u 
solves the Cauchy problem for the wave equation 

u tt - Am = 

with initial data 

u(0,x) = f, u t {0,x)=r 1 Vf. 
Then, as a consequence of the well known decay estimates for solutions to the free 
wave equation (see e.g. [34]), we obtain 

\u(t, X )\ <j(\\f\\ Kl + p/HbJ 

and hence 

(7.3) Ht,x)\<j\\f\\ Ki . 

Here Bf x is the homogeneous Besov space, with norm 

\\v\\ 6ti =J2 2JS \\M^MLi 

jez 

where <j>j now is a homogeneous Paley-Littlewood sequence, i.e., fixed a test function 
il>{r) € Cg such that i[)(r) = 1 for r < 1, ij){r) = for r > 2, we have (f>j(r) = 
i)(2-i+ 2 r) - ijj(2-3 +1 r) for all j e Z. 

The proof of Theorem 1.3 follows the same lines as the proof of Theorem 1.1. 
Consider the Cauchy problem with frequency truncated data 



(7.4) 



iut(t, x) — T>yu(t, x) 
u(0,x) = <Pj(D v )f, 



where (<fij(X))j = o.i,... is the standard Paley-Littlewood partition of the unity defined 
in (5.1). By means of spectral formula, we can represent the solution of (7.4) as 

1 f + °° 

(7.5) U ^^ = 2ViJ ViWe %Rv(\)]fd\. 
Using the identity 

(7.6) R V (X) = R V - RvVR v (I + VR v )-\ 
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which is valid thanks to Corollary 3.7, we can split the integrals in (7.5) into two 
terms, the first one containing the contribution of the free resolvent Rx> and the 
second one containing the contribution of the operator Rx>VRx>{I + VR-d)^ 1 . The 
first term 



A:=^~J ipj{\)e iXt $s [Rt>(\)] f dX 
was estimated above (see (7.3)); it remains to estimate the term 

(7.7) b = -—J_ ^(A)e* At 3 [Q(A)] fdX, 

where 

Q(A) := R V (X)VR V (X)(I + VR v {X)y\ 
After an integration by parts, we obtain 



i p j (X)e iXt —l(Q(X))fdX + J ip' j (X)e iXt T[Q(X)]f dX 
i shows that 

RlVRvih + VRv)- 1 + R v VRl(h + VRv)' 1 



an explicit computation shows that 
dQ 
dX 

+R v VR v (h + VRv^VRKh + VRv)' 1 . 

Now we can apply Lemma 4.2: under assumption (1.21), estimates (4.21), (4.22) 
and (4.23) are satisfied, and the Lemma gives 

(7-9) HQ(A)/|U« <C{Xf\\{ x fl^f\\ L ,, 



(7.10) 



dX 



QWf 



<C(Xf\\( x f/ 2+ cf\\ L2 , 



for some C > 0. Using (7.9) and (7.10) in (7.8) we arrive at the estimate 



\B\<j\\f\\ Ll 



+ 00 



((A) 3 |^(A)| + (A)>;-(A)|) dA 



Recalling that <fijW = 00(2 J 'A), after a change of variables 2 3 A = /i we easily 
obtain 



(7.11) 



\B\< j&Mx 



,3/2+ £/ | 



L 2 • 



From this point on, we can proceed as in the proof of Theorem 1.1 and complete 
the proof of Theorem 1.3. 

Appendix A. Lorentz spaces 

For the convenience of the reader, we recall here the definitions and the main 
properties of the Lorentz spaces L p,q , in view of the applications needed in the proof 
of our results. 

For any measurable function / : K" — > C and any s > we define the upper-level 
£■/ as the set 

E{ : {,• : \f(x)\ > s}. 
The non-increasing rearrangement of / is then the function 

f*(t) := inf{ S > : \E{\ <t}, t G (0,+oo). 
It is also useful to consider the average of /* defined by 

r(t) = \ ff*(r)dr. 
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The standard definition of the Lorentz spaces is the following: 

Definition A.l. For any 1 < p < oo and 1 < q < oo we define the quasinorm 
||/||lp,<j as follows: 

When p ^ 1, if we replace /* with /** in the above definitions we obtain an 
equivalent quasinorm which is actually a norm (see [3], [9]). The Lorentz space L pq 
is defined by 

(A.2) V* = {/ : \\f\\ LP , < oo}. 

Moreover we define 

L 1 ' 1 :— L 1 , L oo,oo = L oo_ 

The spaces L°°' q for 1 < q < oo are usually left undefined (although L 00 ' 1 is defined 
in [9] as the closure of L°° compactly supported functions in the L°° norm). 

With the above definitions, one obtains the elementary properties 

l p,p = LPj 1 < p < oo; 

L M1 CL», l<p<oo, l<qi<q 2 <oo 
(with continuous embedding). When the second index is oo wc obtain the weak 
Lebesge spaces (Marcinkiewicz spaces): 

L P,co = L P^ l<p<oo. 

Moreover, the Lorentz spaces can be obtained by an equivalent construction using 
real interpolation: 

LP , q = {L p» iL pi )f)qi p -i = (i _ 9 ) P - 1 + Bp- X x 

provided 

Po<P\, po<q<oo, < 8 < 1 . 
An alternative characterization of the Lorentz norm can be given using the so- 
called atomic decomposition: 

Lemma A.l. Let f : K™ — > C be a measurable function and let 1 < p < oo, 1 < 
q < oo; then f € L p q if and only if there exist a sequnce of sets (Ej)j e z and a 
sequence of numbers a — (aj)j e z such that \Ej\ = 0(2 : >), a e l q and the following 
estimate 

(A.3) \f(x)\<Cj2aj2- j/p X Ej (x) 

holds, for some C > 0. 

ft is possible to see that the best constant C in (A.3) is equivalent to the Lorentz 
norm of the function /. 

The most useful properties of Lorentz spaces are the Holder and Young inequal- 
ities, which extend the classical ones for Lebesgue spaces. These were originally 
proved by O'Neill in [27]. Wc collect them in the following theorems: 

Theorem A.2 (Holder inequality). Let f e L Pl - q \ g e L P2 - q2 . The following 
estimates hold: 

• ifpi,P2,P e]l,oo[, qi,q 2 ,q & [l,oo], then 

(A.4) \\fg\\ L ^ <C\\f\\ LP1 , n \\g\\ LP 2„ 2 , 1 > p^ 1 = P~\q^ + q 2 1 > 

• ifp\,P2 G [l,oo[, q 1 ,q 2 G [l,oo], then 

(A.5) ||/ff|Ui<C||/||LPi.«iN|iP 2 .«, Pl 1 +to 1 = l, q^+q^^l- 
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We remark that the above statement does not cover the trivial inequality 

(A.6) \\f9\\ L ™ < ||/||i~N|ip„ 

which is easily proved to be true for all cases when L p,q is defined. 

Theorem A.3 (Young inequality). Let f E L Pl ' qi , g e L P2 ' q2 . Then the following 
estimates hold: 

• ifpi,P2,p e]l,oo[, qi,q 2 ,q & [l,oo], then 

(A.7) ||/*ff||LP.«<C , ||/||iPi.«iN|i«.«, Pi'+Pz 1 = l+p-\ gf'+C 1 > 

• ifpi,P2 e]l,oo[, qi,q 2 S [l,oo], then 

(A.8) ||/*ff||Lco<q|/||iPi.,iy|iP 2 .«, p^+p^ = l, q^+q^^l. 

As before, we remark that the above statement does not cover the inequality 
(A.9) \\f*g\\is>.o < C\\f\\ L oo\\g\\ LP „ 

which is easily seen to be true in all cases when L pq is defined (e.g., by real inter- 
polation). 

We conclude this section by studying the weight functions wp(x) — \x\{\ log |a;|| + 
with (3 > 1 which plays a crucial role in our results; in the following proposition 
we determine precisely to which Lorentz the powers w^ s belong. 

Proposition A. 4. For any s > 0, q E [l,oo] we have w^ s E L n / S ^ j provided 
(3 > l/sq. 

Proof. We will use the equivalent Lorentz norm (A.3). For any jgZ consider the 
ball Bi := B 2j/n = {x : \x\ < 2^™} and the rings Ej := LV+ 1 \ B j ; it is clear that 
\Ej\ = C„2 J , where C n depends only on the dimension n. Then, for all x e W 1 we 
have the estimate 

I V(*)l - I g j^fto^ii+l^ ^^l ^ g^'l l0g2 + V-^-^XEM- 

The proof is concluded by the remark that the sequence aj = log2 + l)~^ s is 
in l q if and only if (3 > l/sq. □ 
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